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COHERENT  LASER  RADAR  SYSTEM  THEORY 


Abstract 

Coherent  laser  radars  for  tactical  sensor  applications  are  under 
development  at  a  number  of  laboratories,  based  on  the  mid-infrared 
technology  of  C02  lasers  and  HgCdTe  photodetectors.  Under  U.S.  Army 
Research  Office  Contract  DAAG29-84-K-0095,  a  program  of  research  was 
pursued  to  advance  the  system  theory  of  such  radars,  and  to  corroborate 
these  advances  through  experiments  performed  using  the  test  bed  coherent 
laser  radars  of  the  MIT  Lincoln  Laboratory  Opto-Radar  Systems  Group. 
Toward  those  ends,  fundamental  results  were  derived  for  the  transverse 
and  longitudinal  correlation  scales  of  speckle  targets  observed  via 
heterodyne  detection,  and  pixel-level  statistics  were  derived  and 
experimentally  verified  for  2-D  pulsed  imager  radars  that  use  peak- 
detection  pre-processors.  In  addition,  a  multipixel  multidimensional 
target  detection  theory  was  established  whose  quasi -optimal  processors 
coincide  with  some  ad-hoc  designs  already  in  use,  and  whose  perfor¬ 
mance  analysis  provides  unprecedented  insights  into  the  tradeoffs 
between  radar  system  parameters  and  target-detection  capability. 

Work  was  also  begun  on  the  theory  of  unconventional  laser  radar 
imagers,  e.g.,  synthetic  aperture  systems,  and  on  target-tracking 
theory  for  extended  speckle  objects. 


I.  Research  Summary 

The  development  of  laser  technology  offers  new  alternatives  for  the 
problems  of  target  detection  and  imaging.  Indeed,  coherent  laser  radars  based 
on  the  mid-infrared  technology  of  CC^  lasers  and  HgCdTe  photodetectors  are 
under  development  at  a  number  of  laboratories  [1]  -[5].  The  performance  of 
such  systems  is  strongly  affected  by  the  speckle  patterns  that  are  produced 
by  target  roughness  on  wavelength  scales  [6],  [7].  This  document  is  the  final 
report  on  a  research  program  to  develop  a  quantitative  system  theory  for  such 
radars  through  a  combination  of  analysis  and  experiment.  The  central  issues 
for  this  program  were  the  impact  of  laser  speckle  on  fundamental  pixel  sta¬ 
tistics  [8]  -[10]  and  on  the  design  and  performance  of  multipixel  target 
detection  processors  [11]  -[13].  In  both  cases,  emphasis  was  placed  on  multi¬ 
dimensional,  e.g.,  range  and  intensity,  measurements.  Moreover,  the  funda¬ 
mental  pixel  statistics,  which  served  as  the  foundation  for  the  detection 
analysis,  were  experimentally  verified  [9],  [10]  under  a  collaboration 
arrangement  with  the  Opto-Radar  Systems  Group  of  the  MIT  Lincoln  Laboratory 
using  one  of  their  test  bed  COg  laser  radars  [1].  Finally,  preliminary 
analyses  were  begun  in  the  areas  of  unconventional  laser  radar  imaging  [14] 
and  laser  radar  tracking  theory  [15],  and  laser  reflectometer  measurements 
[16]  were  used  to  support  theory  from  [17].  In  what  follows,  we  shall  sum¬ 
marize  the  principal  results  that  were  obtained  in  the  preceding  problem 
areas. 

Speckle  Statistics  [8]: 

In  order  to  understand  the  impact  of  speckle  fluctuations  on  the  full 


panoply  of  coherent  laser  radar  measurements,  we  derived  the  transverse  and 
longitudinal  degrees  of  coherence  for  speckle  targets  observed  via  heterodyne 
detection.  This  work  elucidated  hitherto  unidentified  interactions  between 


various  measurement-configuration  parameters  that  affect  speckle-target  corre¬ 
lation  scales,  and  hence  laser  radar  performance. 

Peak -Detect ion  Pixel  Statistics  [9],  [10]: 

In  2-D  pulsed  imager  and  2-D  Doppler  imager  radars,  the  intermediate 
frequency  return  signals  are  generally  filtered,  envelope  detected,  thresholded, 
and  peak  detected  in  a  pre-processor  subsystem  [1],  [18].  We  have  analyzed  and 
experimentally  verified  the  resulting  pixel  statistics  produced  with  such 
systems  by  speckle  targets,  and  have  quantified  the  associated  dropout  and 
anomaly  effects  in  range  and  Doppler  measurements. 

Multipixel  Multidimensional  Detection  Theory  [11]  -[13]: 

We  have  addressed  the  somewhat  idealized  problem  in  which  a  2-D  pulsed 
imager  laser  radar  is  used  to  detect  the  presence  of  an  extended  statistical ly- 
uniform  speckle  target  embedded  in  a  statistical ly-uniform  extended  speckle 
background  when  the  target  location  and  target  contrast  are  unknown.  Quasi  - 
optimum  intensity -only ,  range-only,  and  joint  range- intensity  processors  were 
derived,  and  receiver  operating  characteristics  were  computed  for  the  intensity- 
only  and  range-only  cases.  This  work  is  important  because  it  builds  from  the 
correct  pixel  statistics  found  in  [9],  [10],  and  because  the  intensity-only  and 
range-only  processors  coincide  with  ad  hoc  approaches  already  in  use.  The 
work's  greatest  significance,  however,  lies  in  its  quantitative  performance 
predictions,  which  permit  assessing  tradeoffs  between  radar  system  parameters, 
e.g.,  spatial  resolution,  range  resolution,  etc.,  and  target  detection 
performance. 

Target  Reflectivity  Measurements  [16]: 

We  performed  a  series  of  reflectivity  measurements  on  a  variety  of  cal¬ 


ibration  plates  and  spheres  using  an  Incoherent  10.6  ym  wavelength  reflec- 


tometer.  Reflectometer  data  for  the  plates  were  found  to  be  in  close  agree¬ 
ment  with  measurements  collected  with  the  MIT  Lincoln  Laboratory  2-D  pulsed 
imager  laser  radar  test  bed,  as  expected  from  theory  [17]. 


Unconventional  Laser  Radar  Imaging  [14]: 

We  have  been  developing  system  theory  results  for  laser  radar  versions 
of  1-D  and  2-D  synthetic  aperture  radars  (SARs)  and  range-Doppler  (RD)  imagers. 
In  both  cases  our  focus  has  been  to  understand  the  combined  effects  of  target 
speckle  and  local-oscillator  shot  noise  on  system  performance.  The  effects 
of  atmospheric  turbulence  and  laser  frequency  instability  are  also  being 
treated . 

Laser  Radar  Tracking  Theory  [15]: 

We  have  begun  developing  a  theory  for  laser  radar  tracking  of  extended 
speckle  targets.  Thus  far,  we  have  solved  the  track-while-image  problem  in 
which  an  intensity  centroid  estimate  from  the  nth  image  frame  is  used  as  the 
observation  equation  for  a  Kalman-fi Iter  tracker.  The  updated  target  position 
estimate  obtained  from  this  tracker  is  then  used  to  set  the  radar's  optical 
axis  for  the  (n  +  l)st  image  frame. 
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Multipixel,  multidimensional  laser  radar  system  performance 
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United  States  Air  Force  Academy,  Colorado  80840 

Jeffrey  H.  Shapiro 

Department  of  Electrical  Engineering  and  Computer  Science 
Massachusetts  Institute  of  Technology,  Cambridge,  Massachusetts  02139 


Abstract 

The  superb  angular,  range,  and  Doppler  resolutions  of  coherent  laser  radars  have  led 
developers  to  design  imaging  radars  in  multiple  measurement  dimensions.  Designing 
processors  to  detect  targets  in  the  images  generally  proceeds  in  an  ad  hoc  fashion  and  it 
is  difficult  to  predict  the  performance  of  the  resulting  processors.  This  paper  proposes 
simplified  statistical  models  for  the  target,  radar,  and  signals  then  uses  classical 
detection  theory  to  derive  quasi-optimal  processors  which  take  advantage  of  the  aultipixel, 
multidimensional  nature  of  the  image.  The  target  model  is  of  a  radar  looking  down  at  a 
vertical  target  against  a  uniform,  sloping  background.  The  paper  also  presents  the 
receiver  operating  characteristics  (ROCs)  for  the  resulting  generalized  likelihood  ratio 
test  (GLRT)  processors.  The  receivers  may  use  any  combination  of  intensity,  range,  and 
Doppler  measurements.  The  target  reflectivity,  range,  and  angular  location  are  unknown  and 
the  background  reflectivity  is  also  unknown.  The  forms  of  the  quasi-optimal  receivers 
provide  analytical  confirmation  of  the  principles  used  in  many  ad  hoc  processors.  The  ROCs 
not  only  give  bounds  on  the  performance  of  any  ad  hoc  processors  and  prove  the  range-only 
processors  are  usually  superior  to  the  intensity-only  processors,  but  go  on  to  predict  how 
much  better  and  under  what  conditions.  The  ROCs  also  predict  how  performance  changes  as  a 
function  of  resolution  in  one  or  several  measurement  dimensions. 

Introduction 

The  advent  of  laser  sources  with  high  stability,  spectral  purity,  and  sufficient  power 
has  allowed .system  designers  to  translate  much  of  microwave  radar  theory  to  the  optical 
regime.  '  Because  the  spatial  resolution  obtainable  at  optical  and  infrared  (IR) 
wavelengths  is  on  the  order  of  microradians,  most  researchers  have  opted  for  building  and 
analysing  systems  which  perform  a  raster  scan  of  the  target  and  build  an  image  of  the 
target  much  like  a  television  image.  These  radars  are  capable  of  building  intensity, 
range,  or  Doppler  images  or  using  any  combination  of  these  measurement  dimens ions.*' 3 
Although  there  has  bean-much  work  on  analyzing  the  statistics  of  the  target  returns  on  a 
single  pixel  basis,  '  most  of  the  work  on  processing  the  multipixel  images  has  rested 

on  ad  hoc  processors  or  parallel  results  froa.applications  in  the  fields  of  robotics, 
machine  vision,  and  artificial  intelligence.  ' 11  Although  these  approaches  to  image 
processing  have  produced  useful  processors,  they  ignore  the  underlying  statistical  nature 
of  the  image  and  it  is  very  difficult  to  predict  the  processor's  performance  or  how  the 
performance  might  change  as  various  system  parameters  change.  This  work  extends.models  and 
single  pixel  probability  density  functions  (pdfs)  introduced  in  earlier  studies1^  to 
construct  pdfs  for  the  aultipixel,  multidimensional  image  data.  These  density  functions 
(and  the  results  which  follow)  are  suitable  for  radars  which  measure  target  and  background 
reflected  intensity,  range,  velocity  (Doppler  shift),  or  any  combination  of  these 
parameters.  9  (This  paper,  however,  deals  only  with  intensity  and  range.)  The  models  are 
slightly  simplified,  but  still  allow  for  unknown  target  angular  location,  range,  and 
reflectivity  and  unknown  background  reflectivity.  Generalized  likelihood  ratio  test  (GLRT) 
processors  and  their  performance  measures,  including  receiver  operating  characteristics 
(ROCs),  are  derived  for  the  binary  hypothesis  testing  problem.  From  the  ROCs  it  is 
possible  to  predict  the  impact  on  system  performance  from  changing  the  various  system 
psrsmeters  like  radar  optics  aperture,  radar  power,  range  resolution,  and  Doppler 
resolution,  to  name  only  a  few.  These  ROCs  define  the  fundamental  limits  to  the 
performance  of  any  processor  imposed  by  the  statistics  of  the  signals.  They  are 
particularly  useful  as  benchmarks  for  comparing  real  processors.  The  GLRT  processors 
provide  snalytical  confirmation  of  some  of  the  processing  principles  used  in  many  ad  hoc 
approaches. 
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|lim  Iwotkiti  fating 

MU  work  U  eeae«(M<  with  the  problem  of  deciding  whether  or  not  thorn  i«  •  target  in 
•  |if«a  *»1um  of  09000 .  Mis  io  tho  binary  hypothesis  tooting  probloo  t  tbo  radar  oakoa 
woaooraooota  (receives  target  returns)  over  a  volume  of  space  and  uses  the  aeeaureaenta  to 
decide  wbotbor  there  is  a  target  in  the  volume  (bypotnesia  *,  is  true)  or  there  ia  no 
target  ia  tbo  voluae  (hypothesis  sQ  is  true) .  Pigure  1  pictures  the  problem  for  a  voluae 
defined  by  tbo  angular  uncertainty  a  and  range  uncertainty  l  .  The  radar  aakes 
aaaaureaents  ia  the  voluae  by  raster  scanning  the  radar  bean  pattern  acrosa  the  target 
region  dofiaod  by  tbs  angular  uncertainty  a  .  The  radar  oust  discriainate  between  target 
aad  background  while  it  docs  not  know  the  target's  range,  angular  location,  or 
reflectivity. 


Me  laser  radar  nodal  will  be  of  a  connect,  nonostatic,  coherent  laser  radar.  The  radar 
boon  pattern  is  raster  scanned  across  an  uncertainty  region  a  which  nay  contain  a  target. 
Mo  laser  traaaaits  a  series  of  pulses  so  the  received  signal,  after  heterodyning,  is  a 
series  of  either  target  or  background  returns  (depending  on  which  is  illuainated)  plus 
local  oscillator  (LO)  shot  noise.  Sack  target  or  background  return  foras  one  picture 
eleaMnt  (pisel)  of  the  resulting  raster  scanned  iaage.  No  assune  pixels  are  essentially 
non-overlapping  and,  hence,  independent,  be  will  consider  only  purely  speckle  reflectors 
in  this  work.  For  a  ranging  radar,  tbs  pulses  will  be  assuned  to  be  short  duration, 
transfora  linitod  wave for as  which  do  not  resolve  any  range  variations  within  a  pixel.  We 
follow  the  heterodyne  detector  with  a  filter  whose  iapulse  response  ia  notched  to  the 
trsnsaitted  wave fora.  Me  next  step  is  to  square  and  meek  detect  the  filter  output.  The 
peak  detector  outputs  are  then  two  randea  variables,  I,  the  intensity  of  the  peak,  and  t, 
the  oaxiaua  likelihood  (ML)  eatiaate  of  the  reflector  range.  These  outputs  are  our 
aaasured  data  and  tbeir  pdfs  are  known.  Mis  pre-processor  structure  is.shown  in  Figure  2. 
Mis  structure  follows  that  used  in  aany  systens  currently  under  study, **  so  the  results 
are  easily  coapared  with  data  froa  raal  systeas. 

Oacnetty  nodal 

Figure  3  shows  the  aodel  for  the  laser  aad  target  geoaetry.  The  radar  is  above  and 
looking  down  on  the  target  which  is  vertical  to  the  ground.  The  target  angular  subtense  is 
greater  than  a  radar  beaa  width,  so  the  target  is  resolved  in  angle  space.  We  will  call  I 
the  nuaber  of  pixels  on  the  target,  kll  target  pixels  will  occur  at  about  the  sane  range 
since  the  target  is  vertical,  however,  the  background  pixels  will  appear  to  slope  away  froa 
the  radar.  Further,  if  we  know  the  radar's  height  above  the  ground  and  its  pointing  angle, 
we  can  calculate  the  range  to  the  background  if  it  is  reasonably  saooth.  Our  aodel  will 
aasusm  this  range  is  known.  This  aodel  is  psrttcularly  well  suited  for  an  airborne  radar 
looking  for  targets  on  the  ground. 

If  we  know  the  target  site  (angular  extent)  and  shape,  we  can  tile  the  radar  field  of 
regard  with  M  target  shapes  as  in  Figure  4.  we  will  call  one  observation  of  the  radar 
field  of  regard  (MW  pixels)  a  fraae  and  each  target  shape  (M  pixels)  a  subfrane.  we  will 
henceforth  assune  the  targets  align  with  the  subfrane  boundaries  (the  subfrane  contains 
either  an  entire  target  or  no  part  of  a  target)  so  the  subfranes.are  independent  (have  no 
pixels  in  coonon).  Cater  research  has  relaxed  this  assumption. 

single  pixel  density  CwicUqnd 

In  order  to  derive  the  optimal  processors,  we  need  to  know  the  density  functions  for  the 
neasurenents.  Since  the  pixels  are  all  Independent,  we  need  to  find  the  density  function 
for  a  single  pixel.  From  this  density  we  can  generate  nultiplxel  density  functions  easily. 
The  single  pixel  density  functions  for  this  receiver  structure  have  been  published  before 
and  are  repeated  below,  understanding  the  density  and  the  resulting  processors,  however, 
relies  on  understanding  the  operation  of  the  peak  detector  as  well. 

We  can  nodal  the  operation  of  the  peak  detector  by  dividing  the  range  uncertainty  L 
into  0  kins  of  width  equal  to  the  radar  range  resolution.  The  output  of  the  square  law 
envelope  detector  is  approxinately  constant  over  a  bin  tine  and  the  bins  are  approximately 
independent.  For  any  one  bin,  the  intensity  output  is  an  exponential  random  variable  with 
aean  1  if  there  is  no  reflector  at  that  range  or  nean  CWR+1  if  there  is  a  reflector  at  that 
range  (wheae  CNR  is  the  usual  speckle  tsrget  radsr  cart ler-to-noise  ratio  for  the 
reflector17).  Figure  5  shows  the  rsnge  bin  model  with  the  Q  bins  separsted  into  Q.  bins 
known  to  contain  the  background  pixels  snd  Q_  potential  target  range  bins.  The  target 
falls  in  bin  and  the  background  falls  in  Bin  Qw.  The  peak  detector  selects  the  largest 
of  the  Q  intensity  random  variables  and  declares  this  as  the  reflector  intensity  and  the 
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associated  rang*  value  as  tha  rsflsetoc  rang*  for  tha  pixel  in  question.  if  this  procsdura 
salacts  tha  wrong  bin,  wa  say  an  anoaaly  has  occurrad  as  pictured  in  Figure  5. 


Tha  rasulting  aarginal  statistics  pre 
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where  a  ■  (CNR+1)”*,  random  variable  <3  is  the  bin  where  the  peak  occurred,  Qt  is  the  bin 
actually  containing  the  target  or  background  reflector  (where  i  -  t  or  b) ,  and  p  is  the 
probability  of  an  anoaaly  given  by:  * 

PA  ■  *  *  a[  lo9(Q>  -  1/2Q  +  0.577  ]  (3) 

where  TO  is  tha  gasusa  function.  The  approximation  is  valid  for  large  CNR  values.** 

These  results  era  specifically  for  the  ranging  radar.  However,  it  is  possible  to  show 
an  exact  duality  between  all  the  ranging  radar  results  and  those  for  a  Doppler  radar.  3  It 
is  also  possible  to  demonstrate  an  exact. duality  with  a  range  and  Doppler  radar  if  it  has 
an  ambiguity  function  which  is  unimodal.  Henceforth,  all  results  are  for  the  ranging 
radar  with  tha  appropriate  analogies  to  Doppler  systems  understood. 

Binary  detection  receivers 

With  pixal  statistics  in  hand,  we  can  derive  the  optimal  receiver  for  choosing  between 
the  two  hypotheses  H.  and  H..  We  will  use  the  Neyman-Pearson  criterion  which  constrains 
F««  tha  probability  of  falsa  alarm,  to.be  less  than  or  equal  to  a  specified  value  and 
minimises  PM,  the  probability  of  miss.  *  The  result  is  a  likelihood  ratio  test1*: 


A(f) (r) 
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where  the  superscript  (f)  emphasises  the  data  vector  r  and  the  likelihood  ratio  a(£>()  are 
the  measurement  data  and  tha  likelihood  function  for  the  entire  frame.  Since  all  pixels 
give  independent  measurements  (minimal  beam  overlap)  and  tha  target  is  entirely  in  one 
subframe,  tha  frame  density  function  is  simply  the  product  of  the  subframe  density 
functions  which  are,  in  turn,  the  product  of  the  pixel  density  functions.  The  threshold 
is  chosen  to  meet  the  Pr  constraint  with  equality. 

gafrrerc  FfgWtWt 

Tha  density  functions  introduced  in  the  last  section  depend  on  unknown,  non-random 
parameters  like  tha  target  or  background  range  bin,  Qt .  we  eliminate  the  unknown 
parameters  by  using  the  generalised  likelihood  ratio  test  (GLRT) : ** 
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where  A  is  tha  unknown  parameter  vector.  For  our  problems,  the  unknown  parameters  are  m  , 
tha  actual  subframe  containing  tha  target i  Qt,  tha  actual  target  range  bin  number;  and  b8th 
target  and  background  CNRsi  CNR^  and  CMRg,  respectively. 

Introducing  tha  GLRT  hare  is  a  crucial  step  in  the  development.  The  unknown  parameters 
involve  tha  aultipixel  nature  of  the  target  and  allow  us  to  extend  single  pixel  statistical 
analyses  to  tha  aultipixel  case. 


We  can  make  additional  simplifications  to  the  frame  likelihood  ratio.  Since  subframes 
do  not  overlap  and  tha  target  is  always  aligned  with  a  subframe  (by  assumption)  we  can 
separate  the  frame  density  function  into  subframe  density  functions  which  only  differ, 
under  hypotheses  H0  and  H.,  for  the  one  subfrang  a  .  we  can  further  simplify  the  density 
functions  because  the  unknown  parameter  vector  A  contains  some  elements  which  affect  the 
density  under  hypothesis  Hq  but  not  under  H^  and  vice  versa.  This  fact  allows  us  to  make  a 
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siaplification  foe  largo  valuta  of  N  (tht  nuabor  of  subfraaos  par  fraat) .  Tho  raault  iat 


max  { A 1**  } 
l<ra<M  gm 


£s(s) iH(a)  J  (r‘8) 
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where  tho  (a)  superscript  indicates  a  subfraao  quantity  and  tho  a  subscript  is  still  a  sub- 
frame  indox.  (Hypothesis  H,'  1  indicates  a  target  present_in  this  subfraao  and  hypothesis 
H.1  '  indicates  no  target  in  this  subfraao).  The  vectors  A.  and  A.  are  the  portions  of  the 
paraaeter  vector  A  affecting  the  density  under  hypotheses  Hq  and  Hw  respectively.  The  * 
indicates  an  ML  estiaate  of  the  paraaeter  subvector.  u  1 


Because  of  the  fora  of  Equation  (6)  it  is  easy  to  show  the  fraae  level  statis 
and  Pg  depend  on  the  subfraao  statistics  Pp'  '  and  PM(  '  in  a  siaple  fashion: 
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where  the  approx iaat ions  are  valid  for  the  usual  PF(>><<1  case. 

These  equations  have  iaportant  physical  interpretations.  The  false  alara  probability 
rises  linearly  as  the  nuaber  of  subfraaos  (the  angular  search  area)  increases.  The  aiss 
probability  is  approxiaately  independent  of  the  search  area.  A  detection  occurs  if  any 
subfraao  statistic  clears  the  threshold,  even  if  it  is  not  the  correct  subfraao.  The 
probability  of  detection  on  the  wrong  subfraao  is  quite  saall,  however,  since  this  is 
basically  a  subfraao  false  alara  with  probability  Pp*  . 

Intensitvonlv  processors 

First  consider  the  processor  which  uses  only  the  aeasured,  peak  detected  intensity  for 
each  pixel.  Because  of  the  peak  detector  the  pdf  for  i  is  coaplicated  and  it  is  difficult 
to  derive  an  exact  optiaal  processor.  To  derive  the  processor,  we  used  a  Central  Liait 
Theorea  approximation  to  the  density.  This  gives  an  eainently  reasonable  processor  as  we 
shall  shortly  see.  To  analyse  the  processor  perforaance,  there  are  better  approx iaat ions 
to  the  tails  of  the  density  function  (based  on  aodified  Chernoff  bounds)1*  which  give  aore 
accurate  results  than  the  Central  Liait  Theorea. 

Applying  the  Central  Liait  Theorea  approx iaat ion,  perforaing  the  required  algebra,  and 
siaplifying  the  expression,  ^  the  final  log  likelihood  ratio  for  the  quasi-optiaal 
intensity-only  processor  is: 

log  A***  »  max  {log  a!** } 

^  l<m<M  9,0 


max  { |  V  I. 
l<m<M  n«l 


Nubl> 


where  I„  is  the  aeasured  intensity  for  the  n-th  pixel  in  the  a-th  subfraae,  and  uh 
ML  estilate  of  the  background  aean  intensity  under  hypothesis  H0:  ° 


is  the 
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The  processor  in  Equation  (10)  ha/  s  reassuring  fora.  It  says  we  aake  an  estiaate  of 
the  average  background  intensity  and*coapare  it  to  the  average  intensity  of  each  subfraae. 
The  subfraae  whose  intensity  is  aost 'different  froa  the  background  is  declared  the  target 
subfraae  if  the  difference  is  greater  than  the  threshold.  If  the  difference  is  less  than 
the  threshold,  we  declare  flQ:  no  target  present.  In  other  words,  the  processor  searches 
for  target-to-background  intensity  contrast.  This  is  exactly  what  intuition  tells  us  to  do 
and  what  aany  researchers  have  done  with  their  ad  hoc  processors. 


Now  consider  a  processor  which  uses  only  the  range  bin  inforaation,  q  ,  the  aeasured 
bin  nuaber  where  the  peak  intensity  occurred,  for  this  processor  we  use*ehe  single  pixel 
range  statistics  given  earlier,  introduce  the  aultipixel  statistics  by  using  pixel 
independence,  and  aake  GLUT  processors  by  separating  the  unknown  paraaeters  under 
hypotheses  Hq  and  H^.  The  approxiaate  generalised  likelihood  ratio  which  results  is: 


max  {  A'®' 
l<m<M 


}  »  max  {  j  -  k  } 
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where : 
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Here  Q.  is  the  (known)  range  bin  nuaber  for  the  background  in  pixel  a,  nj  q  is  the 
peak  detector  output  bin  nuaber  for  pixel  a,  n;  and  fl.  is  the  presumed  target  range  bin 
nuaber  for  the  a-th  subfraae.  In  words,  k  is  the  number  of  tiaes  the  peak  detector  found 
the  peak  intensity  occurred  at  the  cor rect^iac kg round  range  in  the  a-th  subfraae.  Randoa 
variable  j  is  the  nuaber  of  tiaes  the  peak  detector  found  the  intensity  peak  at  the 
presuaed  target  range  bin  in  the  a-th  subfraae.  The  presuaed  target  range  bin,  for  a 
particular  subfraae,  is  the  potential  target  range  bin  which  the  peak  detector  chooses  aore 
often  than  any  other  potential  target  range  bin  for  that  subfraae. 

Physically,  Equation  (12)  says  the  processor  deteraines  whether  the  presuaed  target 
range  bin  or  the  known  background  range  bin  was  selected  aore  often  by  the  peak  detector 
for  each  subfraae.  The  aore  often  the  presuaed  target  bin  is  selected,  relative  to  the 
background  bins,  the  larger  the  statistic  A**’.  If  the  aaxiaua  of  these  statistics  over 
all  H  subfraaes  exceeds  the  threshold,  the  processor  declares  H,:  target  present.  In 
other  words,  the  processor  looks  for  the  range  aeasureaents  to  aggregate  or  cluap  in  either 
the  known  background  range  bins  (if  no  target  is  present)  or  one  of  the  Q_  potential 
target  range  bins  (if  a  target  is  present).  Colloquially  put,  the  processor  searches  for 
range  contrast.  In  particular,  owing  to  the  nature  of  the  Figure  3  geoaetry,  this  range 
contrast  can  also  be  called  verticality. 

Receiver  perfotaance  result! 

It  is  possible,  at  least  in  theory,  to  find  perforaance  aeasures  for  the  processors  just 
derived  by  integrating.the  pdfs  over  the  proper  ranges,  set  by  the  detection  threshold,  to 
get  F_  and  PM  values.4  In  practice  it  is  difficult  or  impossible  to  carry  out  these 
integrationsaxcept  numerically  or  with  approxiaation  techniques.  In  this  work  we  carried 
out  the  discrete  range-only  processor  analyses  numerically  and  the  continuous  intensity- 
only  processors  analyses  using  approximations  derived  froa  the  Chernoff  bounds1*'1  . 


The  processor  equations  contain  aany  paraaeters:  N,  the  nuaber  of  subfraaes  per  fraae; 
M,  the  nuaber  of  pixels  per  subfraae;  Q,  the  nuaber  of  range  bins  per  pixel;  CNR. ,  the 
target  CNR;  CNR.,  the  CNR  of  the  background  pixels;  and  the  contrast  ratio, .L  =  CNR./CNR.. 
Some  of  these  paraaeters  depend  on  each  other  through  the  radar  equation,  '17  so  weccannot 
change  then  without  accounting  for  the  effect  on  other  paraaeters  (aost  notably,  the  CNRs) . 
It  is  important  to  understand  these  relationships  because  we  want  to  coapare  radars  with 
different  resolution  capabilities  iaaging  the  saae  target  and  background  environaent.  If 
we  do  not  account  for  these  paraaeter  Interdependences,  the  coaparisons  are  not  valid. 

For  siaplicity,  we  examine  the  radar  equation  with  a,  the  atmospheric  extinction,  ■  0, 
and  optical  efficiency,  e  ■  1,  so  the  radar  equation  reduces  to1*'17: 


!xis 
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t  or  b) 


where  n  is  the  detector  quantua  efficiency,  hv  la  the  optical  photon  energy,  o.  is  the 
speckle  reflector  diffuse  reflectivity'  (which  Is  constant) ,  P_  is  the  transmitter  peak 
potter ,  B  is  the  aatched  filter  bandwidth,  A_  is  the  radar  optica  area,  and  L,  is  the 
reflector  ranqe.  The  subscript  i  can  be  either  t,  for  the  target,  or  b,  for  the 
background . 

If  we  change  A_,  the  radar  aperture  area,  we  will  change  N,  the  nuaber  of  pixels  on  the 
target.  Both  P_,  the  radar  peak  transmitted  power,  and  B,  the  radar  bandwidth,  are  related 
to  Q,  the  nuabef  of  range  bins  per  pixel,  and  also  whether  the  laser  itself  is  peak  power 
Halted  or  average  power  limited. 

If  we  define  a  CNR  at  aaxiaua  radar  aperture,  and  aaxlaua  bandwidth,  b*“*x* : 
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we  can  find  equations  for  CNR,  at  various  N  and  Q  values  for  peak  (P_)  and  average  (P.„) 
power  Halted  lasers:  AV 
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The  value  Q*“*x*  exists  because  if  Q 


se  if  Q.aets  too  large,  we  would  violate  the  asauaptton  of 
a  ragae  unresolved  target.  The  value  N'  '  siaply  repreaents  the  aaxiaua  aperture  site, 
Ajj*  ,  for  any  given  system  design. 

In  all  the  calculations  that  follow,  we  set  N^a*x*  •  40  pixels,  Q*“x*  •  10,000  bins, 
and  H  •  1000.  (Changing  radar  parameters  does  not  affect  M.) 


[ntensity-on] 


cocesaor 


First  we  will  exaaine  the  intensity-only  processor.  Equation  (10),  dependence  on  target 
CNR  then  look  at  the  dependence  on  Q  and  N.  In  each  case  we  plot  PH  as  a  function  of  the 
one  variable  of  interest  and  keep  all  the  others  constant.  Figure  V  plots  PM  versus  CNR. 
for  N  ■  10  (part  a)  and  N  *  20  (part  b) .  In  each  case  a  plot  is  shown  for  contrast  ratios 
of  c  ■  +5  dB  and  >10  dB  and  for  two  P.  values,  P.  ■  10~J  and  P_  ■  10~  .  M  and  Q  are 
constant  at  N  •  1000  and  Q  ■  10.  r  1 

There  are  two  important  aspects  to  these  plots  which  bear  explaining.  First,  as  CNR 
gets  large,  the  PM  approaches  a  non-zero  asyaptote.  This  occurs  because  we  approach  the 
speckle  Halted  performance  regiae.  At  high  CNR.  values  there  is  essentially  no  LO  shot 
noise  to  contend  with  and  perforaance  is  Halted  by  the  speckle  induced  fluctuations  in 
target  and  background  intensities.  Second,  for  large  positive  contrasts,  the  perforaance 
iaproves.  For  negative  contrasts,  even  large  ones,  the  perforaance  is  very  poor.  For  the 
no  contrast  case  it  is  easy  to  understand  why  the  processor  can  do  no  better  than  to  guess: 
it  is  looking  for  contrast  and  there  is  none.  The  poorer  performance  at  large  negative 
contrasts  than  at  large  positive  contrasts  is  a  result  of  the  asyaaetry  of  the  density 
functions  for  the  intensities  under  the  two  hypotheses.  For  negative  contrasts,  the  target 
aean  intensity  is  lower  than  the  background  mean  intensity.  However,  the  speckle  induced 
fluctuations  in  the  background  aeasureaents  often  give  Intensities  less  than  the  background 
aean  and  closer  to  the  target  aean  intensity.  This  makes  the  two  virtually  impossible  to 
distinguish  at  realistically  low  Pp  values. 

Figure  7a  shown  PM  as  a  function  of  Q,  the  nuaber  of  range  bins,  for  a  constant  average 
power  laser  while  Figure  7b  shows  the  saae  inforaatlon  for  a  constant  peak  power  laser.  In 
each  case  plots  are  for  two'values  of  N,  10  and  20.pixels,  and  two  values  of  CNR.T  ,  14 
and  20  dB,  with  fixed  values  of  M  ■  1000,  P-  ■  10  ,  and  ;  •  >10  dB.  The  figures 
incorporate  the  CNR  corrections  for  chang inf  N  and  Q  values. 

For  the  constant  peak  laser  power  model,  the  perforaance  falls  at  high  Q  values  since 
the  CNR  depends  on  Q.  As  Q  increases,  the  CNR  drops  until  it  moves  away  froa  the  speckle 
Halted  perforaance  regiae.  The  weak  Q  dependence  in  the  constant  average  laser  power 


aodel  perforaance  is  due  to  ths  weak  dspsndsncs  of  anoaaly  probability  P. ,  on  Q  shorn  in 
aquation  (3).  As  Q  increases,  ths  probability  of  anoaaly  also  increases'^  lowly  which 
eausas  a  drop  in  parforaanca. 

Piqura  8a  plots  PM  varsus  H.  ths  oiaber  of  pisals  par  subfraae,  for  a  constant  average 
power  laser,  while  Ffguxe  8b  plots  tffb  saaa  quantities  for  a  constant  peak  power  laser. 

The  plo|f  are  for  values  of  Id-  and  20  dR  and  Q  values  of  10  and  1000  bins.  Again, 

F_  •  10”*  and  N  •  1000  subfraaes.  Ths.plots  account  for  the  CMt  variations  with  M  and  Q. 
Proa  the  fora  of  the  Cbernoff  bounds,  we  expect  the  parforaanca  to  iaprove 

exponentially  for  increasing  M  values  and  the  figures  bear  this  out.  Ths  differences  in 

the  two  plots  are  due  to  the  differences  in  perforaance  for  changing  Q  values  between  the 

two  laser  aodels.  These  effects  we  have  just  discussed  above. 

Ranee-only  processor 

Now  we  exaaine  the  range-only  processor  of  Equation  (12).  We  will  exaaine  the  saae 
variables  we  did  for  the  intensity-only  processor.  Figure  9  plots  versus  CMS.  for  M  ■ 

10  (part  a)  and  W  •  20  (part  b) .  Three  contrast  ratios  are  shown  ana  two  P_  values,  while 

0  •  10  and  N  ■  1000  are  constant.  Here  the  perforaance  laproves  as  CMRt  gets  larger  and 

does  not  bottoa  out  at  a  non-sero  asyaptotic  value.  As  CNRfc  gets  larger,  P.  goes  to  xero 

and  there  are  no  anoaalies  to  confuse  the  processor,  so  it  can  perfora  arbitrarily  well. 

We  also  notice  for  a  fixed  CNR*,  perforaance  gets  poorer  as  contrast  increases.  This 
occurs  because  as  c  increases  ror  a  fixed  CNRfc,  CUR.  aust  decrease.  As  CNRw  decreases,  the 
processor  has  sore  difficulty  distinguishing  Che  background  and  perforaance  falls  off. 

Figure  10  plots  PM  versus  0  for  constant  average  laser  power  (part  a)  and  constant  peak 
laser  power  (part  b)7  For  the  constant  average  power  laser,  the  perforaance  Increases  as  Q 
increases  for  a  while  then  starts  to  fall  slowly.  This  is  a  result  of  two  coapetlng 
factors.  As  Q  increases,  the  probability  of  getting  anoaalies  (which  occur  in  randoa  bins) 
to  cluap  together  well  enough  to  aasquerade  as  a  target  falls.  This  iaproves  perforaance. 
However,  as  Q  increases,  P.  also  increases  approxiaately  logar ithatcally  as  shown  in 
Equation  (3).  This  causes*per f oraance  to  fall.  For  the  constant  peak  power  laser,  the 
increase  in  CNR_  as  Q  falls  overwhelas  the  weaker  Q  dependence  in  P..  Perforaance  laproves 
quickly  as  Q  decreases. 


Perforaance  conoat Isons 

There  are  aany  ways  to  present  the  perforaance  data  and  the  nuaber  of  variables  we  have 
available  aakes  it  difficult  to  present  aore  than  a  saall  portion  of  the  data  at  one  tine. 
Figure  11  shows  a  way  to  present  the  data  that  could  be  particularly  useful  to  a  systea 
engineer.  In  this  plot  we  select  specified  values  for  P.  and  PH  at  a  certain  cm.  and  C. 
Then  we  check  to  see  if  the  range-only  or  intensity-only  processor  aeets  the  perforaance 
criteria  for  various  w  and  0  values.  These  plots  then  indicate  the  nlniaua  coablnatlons  of 
angular  and  range  resolutions  required  to  aeet  the  perforaance  specifications. 


Figure  11,  for  the  constant  average  laser  power  aodel,  shows  the  aore  significant  iapact 
of  laser  aperture  sise,  A.,  (reflected  in  the  N  value)  on  perforaance  relative  to  the  range 
resolution  paraaeter  Q.  Generally  the  processor  either  aeets  requlreaents  for  a  given  M  or 
not,  regardless  of  Q.  Figure  11a  is  fox. relatively  high  perforaance  requlreaentst  P_  • 
99.9%  at  P.  -  10”*  and  fairly  low  CNR  values.  At  1C  dB,  we  find  the  point  N  -  20 
pixels,  Q  •  10  bins  where  the  intensity-only  processor  out-perforas  the  range-only 
processor.  At  higher  Q  values,  _>  300,  the  intensity-only  processor  falls,  but  .the  range- 
only  processor  satisfies  the  requlreaents  for  Q  >  30.  Xf  we  increase  the  CNR*  by  only  4 
dl,  we  can  aeet  requlreaents  with  either  processor  as  long  as  N  _>  20. 

In  Figure  lib  we  relax  the  perforaance  requlreaents  soaewhat  to  PD  ■  95%  at  PF  ■  10-3, 
but  we  reduce  the  contrast  to  only  5  <J».  Sere  the  intensity-only  processors  need  at  least 
N  •  30  pixels  to  aeet  the  relaxed  requlreaents  because  they  perfora  so  poorly  at  low 
contrasts.  The  range-only  processors  perfora  well  enough,  however,  that  for  a  4  da 
increase  in  CNRtl  ,  we  can  cut  angular  resolution  in  half,  froa  N  •  20  pixels  to  N  ■  10, 
pixels  and  still  aeet  requlreaents. 


Figure  12  is  a  presentation  Identical  to  Figure  11,  but  now  for  a  laser  of  constant  peak 
power.  Eere  the  processors  both  perfora  better  at  low  Q  values  and  poorer  at  high  Q  values 
because  the  CNR  rises  for  saaller  Q  values,  as  shown  in  Equation  (1C).  Now  it  is  possible 
for  the  range-only  processor  to  aeet  the  perforaance  criteria  at  lower  angular  resolution 
(N  values)  than  before  because  of  the  Increased  CNR  at  low  Q  values. 


Conclusions 

In  this  paper  we  have  seen  how  it  is  possible  to  derive  nearly  optiaal  processors  for 
the  binary  hypothesis  testing  problea  and  use  the  statistics  to  find  the  perforaance  of  the 
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processors.  Hi  found  ths  qussi-optiasl  processors  are  such  like  those  researchers  have 
used  on  an  ad  hoc  basis.  Although  this  paper  presented  only  the  ranging  radar  results, 
exact  dualities  exist1J  which  apply  the  results  directly  to  radars  eeasuring  intensity, 
range,  Doppler,  or  any  combination  of^  these  three.  The  processors  were  derived  froa 
realistic  aodels  assuaing  nuaerous  qrfknown  parameter a:  target  range,  target  angular 
position,  target  reflectivity,  and  bickground  reflectivity.  Since  it  is  possible  to 
predict  the  performance  of  these  processors  analytically,  we  were  able  to  show  how 
perforaance  varies  as  a  function  of  target  sise,  angular  resolution  or  radar  aperture, 
range  resolution,  target-to-background  contrast,  and  overall  CNR.  A  few  case  studies  were 
presented  as  exaaples.  These  results  aake  it  possible  for  a  systea  engineer  to  select 
appropriate  design  criteria  for  a  laser  radar  and  predict  how  changes  in  one  parameter 
interact  with  other  paraaeters.  They  are  also  useful  because  they  deaonstrate  the  limits 
to  systea  performance  dictated  by  the  noise  and  speckle  statistics.  These  results  are 
useful  as  benchaarks  against  which  other  processors,  even  ad  hoc  ones,  can  be  measured. 
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H,  •'  Target  present 
Hq:  Target  absent 
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Figure  2.  Pra-proceasor  structure  for  a 
ranging  radar. 


Figure  1.  The  target  detection  problem. 


=  12  subframes 
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Figure  3-  Laser-target  geometry  model. 


Figure  4.  Target  subframe  model. 
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Torgst  CNR,  dB 


Figure  5.  Range-bin  model  for  peak  Figure  6  (a).  Intensity-only  processor  miss 

detection.  probability,  Pj|,  vs.  target 

carrier-to-noise  ratio,  CNRt, 
for  N*10  pixels  per  subframe, 
M"1000  subframes  per  frame,  and 
Qa10  range  bins.  Upper  two 


Target  CNR,  dB 


Figure  6  (b).  Intensity-only  processor  miss 
probability.  Pm,  vs.  target 
carrier-to-noise  ratio,  CNR,., 
for  N-20  pixels  per  subframe, 
M*1000  subframes  per  frame,  and 
Q*10  range  bins.  Upper  two 
curves  are  for  C"5  dB  target 
contrast ;  lower  two  curves  are 
for  C*10  dB  target  contrast. 
Dashed  curves  are  for  false- 
alarm  probability  Pp-lQ-°; 
solid  curves  are  for  Ppa10"3. 


curves  are  for  c«5  dB  target 
contrast ;  lower  two  curves  are 
for  c«10  dB  target  contrast. 
Dashed  curves  are  for  false- 
alarm  probability  Pp«10-6;  _ 
solid  curves  are  for  Pp-10_'J. 


♦  Rang*  Bins:  0 


Figure  7  (a).  Intensity-only  processor  miss 
probability.  Pm,  vs.  number  of 
range  bins,  3,  for  constant-PAy 
laser  model  with  M*1000  sub- 
frames  per  frame,  Pp«10"°  false- 
alarm  probability,  and  c«10  dB 
target  contrast.  Upper  two 
curves  are  for  N«10  pixels  per 
subframe;  lower  curves  are  for 
N*20  pixels  per  subframe.  Dashed 
curves  are  for  CNR^'*20  dB; 
solid  curves  are  for 
CNR(j?)-l6  dB. 
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7  (b).  Intensity-only  processor  miss  Figure  8  (a), 
probability,  PM,  vs.  number  of 
range  bins,  Q,  for  constant-Pij 
laser  model  with  M-1000  sub- 
frames  per  frame,  Pp-10"6  false- 
alarm  probability,  and  5-10  dB 
target  contrast .  Upper  two 
curves  are  for  N-10  pixels  per 
subframe;  lower  two  curves  are 
for  N-20  pixels  per  subframe. 

Dashed  curves  are  for  CNR'jjP- 
20  solid  curves  are  for 
CNRtpjMl.6  dB. 
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xncensicy-oniy  processor  miss 

Pr°5ability»  PM*  v»-  number  of 
subframe  pixels,  N,  for 

constant-P»v  laser  model  with 
M-1000  subframes  per  frame, 
Pp-10-6  false-alarm  probability, 
and  c-10  dB  target  contrast. 
Arrows  indicate  curves  for  Q-10 
and  Q-1000  range  bins.  Dashed 
curves  are  for  CNRW.20  dB: 
solid  curves  are  fof  CNR^*'- 
16  dB.  t 


♦Pixels:  N  Target  CNR,  dB 


Figure  8  (b) 


Intensity-only  processor  miss 
probability,  p^,  vs.  number  of 
subfraae  pixels,  N,  for 
conscant-Pip  laser  model  with 
M-1000  subframes  per  frame, 
Pp-10"O  false-alarm  probability, 
and  5-10  dB  target  contrast. 
Arrows  indicate  curves  for 
Q-10  and  Q-1000  range  bins. 
Dashed  and  solid  Q-1000  curves 
CNRCgT.20  dB  and 
CNRlg)-l6  dB,  respectively; 
®*1?o^urv®  is  for  both 
CNR^^-20  dB  and  CNR<?)-16  dB. 
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Range-only  processor  miss  proba¬ 
bility,  Pjj,  vs.  target  carrier- 
to-noise  ratio,  CNRt,  for  N-10 
pixels  per  subframe,  M-1000  sub- 
frames  per  frame,  and  Q-10  range 
bins.  Rightmost  two  curves  are 
for  c-10  dB  target  contrast; 
middle  two  curves  are  for  5-5  dB 
target  contrast;  leftmost  two 
curves  are  for  5*0  dB  target 
contrast.  Dashed  curves  are  for 
false-alarm  probability  Pp«10-6. 
solid  curves  are  for  Pp-lu-3, 
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Figure  9  (b).  Range-only  processor  miss  Figure  10(a). 

probability,  Pj|,  vs.  target 
carrler-to-noise  ratio,  CNRt, 
for  N*20  pixels  per  sub f reuse, 

M“1000  subframes  per  frame, 
and  Q"10  range  bins.  Rightmost 
two  curves  are  for  c«10  dB  tar¬ 
get  contrast;  leftmost  curve  Is 
for  c*5  dB  target  contrast. 

Dashed  curves  are  for  false- 
alarm  probability  Pp«10“°; 
solid  curve  is  for  Pp*10“3, 


Range-only  processor  miss  prob¬ 
ability,  P|y|,  vs.  number  of 
range  bins,  Q,  for  constant-PAV 
laser  model  with  M-1000  sub- 
frames  per  frame,  Pp«10“° 
false-alarm  probability,  and 
C*10  dB  target  contrast.  Upper 
two  curves  are  for  N«*10  pixels 
per  subframe;  bottom  curve  is 
for  N“20  pixels  per  subframe. 
Dashed  curve  is  for  CNR^g^- 
20  dB:  solid  curves  are  for 
C.'JRCOj.ifi  dB. 
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Figure  10  (b).  Range-only  processor  miss  prob¬ 
ability,  P m,  vs.  number  of  range 
bins,  Q,  for  constant-P^  laser 
model  with  M"100Q  subframes  per 
frame,  Pp«10"°  false  alarm  prob¬ 
ability,  and  cs10  dB  target  con¬ 
trast.  Leftmost  two  curves  are 
for  N«10  pixels  per  subframe; 
rightmost  curve  is  for  N*20  pixels 
per  subframe.  Dashed  curve  is  for 
CNR'2'*20  4B;  solid  curves  are 
for  CNRtg'-l6  dB. 
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Figure  11.  (a)  Target-detection  performance  trade-offs  for  constant-P^y  laser  model  with 
*■1000  subframes.  High-performance/high-contrast  case:  "R"  ■  range-only 
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Figure  11.  (b)  Target-detection  performance  trade-offs  for  constant-P^y  laser  model  with 
M-1000  subframes.  Low-performance/low-contrast  case:  "R"  -  range-only 
processor  achieves  Pjh<0.05  at  Pp  •  10"’  for  (  »  5  dS  target  contrast; 

"I"  •  Intensity-only  processor  achieves  Pm  <0.05  at  Pp  -  10"3  for  ;  -  5  dB 
target  contrast. 


Pigure  12.  (a)  Target-detection  performance  tradeoffs  for  constant-Pf  laser  model  with 

M  -  1000  subframes.  High-performance/hlgh-contraat  caae :  "R"  •  range-only 

processor  achieves  Pjh<10"3  at  Pj,  •  io~®  for  ;  -  10  dB  target  contrast ; 

"I"  ■  intenalty-oniy  pl-ocessor  achieves  P^<10“'J  at  Pp  *  10”“  for  c  *  10  dB 
target  contrast. 


log(Q) 


log(Q) 


Pigure  12.  (b)  Target-detection  performance  trade-offs  for  constant-P?  laser  model  with 
M  ■  1000  subframes.  Low-performance/low-contrast  case:  "R"  *  range-only 
processor  achieves  P,m  £  0.05  st  Pc*  ■  10*3  for  c  *  5  3B  target  contrast; 

"I"  ■  intensity-only’  processor  achieves  Pj|  <.  0.05  at  Pp  ■  10“3  for 
;  *5  dB  target  contrast. 
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